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Abstract 

We examine the generalized leading-logarithmic approximation (LLA) equations for 
compound states of n-reggeized gluons. It is shown that in multi-color QCD, when 
Nc ^ oo, these equations have a sufficient number of conservation laws to be exactly 
solvable. Holomorphic factorization of the wave functions is used to reduce the corre- 
sponding quantum mechanical problem to the solution of the one- dimensional Heisenberg 
model with the spins being the generators of the Mobius group of conformal transforma- 
tions. 



1. Introduction 

This talk is centered around obtaining the exact solution to a perturbative QCD evo- 
lution equation known as the Bartels-Kwiecinski-Praszalowicz (BKP)-equation^ in the 
limiting case where the number of colors of gluons, Nc, is infinite. One may wonder what 
relevance any equations of perturbative QCD may have in understanding the low energy 
confining properties of the theory. The answer to this question is not well defined. How- 
ever, what is clear is that only in the perturbative regime of QCD, we are able to exactly 
treat gauge and Lorentz invariance. Even then, within this regime one discovers that such 
a task is nontrivial. Thus the first lesson one gains from examination of perturbative QCD 
is experience with nonabelian gauge calculations that can be tested for their correctness. 

That may be a useful reason for those working in low energy QCD to nevertheless 
study the high energy regime as a warm-up exercise. However that is not the primary 
reason for this talk. The general class of equations that we are considering here are 
the only known evolution equations in QCD that exactly respect gauge invariance and 
have a kinematic regime in which they are exactly valid. Although their practical uses 
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are for calculating high energy scattering amplitudes, it is natural to also examine what 
properties of these equations and their solutions correspond to what we believe to be true 
at low energy. 

To start with, let us first introduce the names of the equations to which we are referring 
and give some history on their development. The first evolution equation in this class 
was the Gribov-Lipatov-Altarelli-Parisi (GLAP) equation^ which was derived in the early 
70's. This equation is a predecessor to the main equations we want to discuss here. From 
this group, the first was the Fadin-Kuraev-Lipatov (FKL)-equation^, which was derived 
in 1975. This equation was the initial form of a more contemporary version known as 
the Balitsky- Fadin-Kuraev-Lipatov (BFKL)-equation'^. The FKL-equation was derived 
for massive Yang Mills theory with a massive Higgs particle and for arbitrary SU(N) 
gauge group. This is an equation for the two-to-two scattering amplitude in the Regge 
limit, m^e^/^ ~ s ^ m^ ~ t, where m is the mass of the vector boson, ^/s is the center 
of mass energy and \f—i is the momentum transfer. We note that the Regge limit also 
implies the leading-log-approximation, where g^ ln(s/t) ~ 1 and g"^ <C 1. To obtain the 
equation, working in momentum space using s-channel unitarity along with analyticity, 
the amplitude was computed to eighth-order. From this the general form could be deduced 
into what became the FKL-equation. Only the t=0 solution was obtained in [3]. The 
solution showed that the amplitude violated the Froissart bound. However it should be 
realized that the region where this violation occurs is also beyond the region where the 
FKL-approximations are valid. 

In 1978 it was verified that there are no infrared divergences in QCD for scattering 
of colorless particles at arbitrary t in the BFKL-equation. In particular this held at t=0, 
where the problem is typically most pronounced. The solution for arbitrary t was found 
in 1986*^. A relevant point for the present discussion is that the calculation was done in 
transverse coordinate representation (or impact parameter space). In this representation 
it was recognized that the BFKL-equation had two-dimensional Mobius invariance. 

The shortcoming of the BFKL-equation is that it violates the unitarity bounds. To 
correct this, the suggestive approach is to consider diagrams with an arbitrary num- 
ber of reggeized gluons. The BFKL-equation only accouts for two reggeized gluons. 
The equation with N reggeized gluons was obtained by Bartels and by Kwiecinski and 
Praszalowicz^ . The purpose of this talk is to examine the solutions of this equation for 
Nc ^ oo. What will be achieved here is a relation of this equation in the above limit, to 
exactly solvable models. The end result is a reduction of the problem to a one-dimensional 
lattice model. 

Before turning to the quantitative discussion, let us place into perspective what contact 
this development makes with the problem of confinement. We have believed since the 
early seventies that Yang-Mills theory is plausibly the low-energy limit of an appropriate 
string theory. In the high energy limit, one may then ask if any aspect of QCD's string- 
like nature manifests. There is no known reason from general principles to expect this. 
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Nevertheless, in light of the results we discuss here, we do find a string-like remnant of 
QCD in this limit. 

Examining the issue a little further, we next recall that high energy processes in 
fact have an intrinsic dependence on the low energy properties of QCD. We have known 
since the 60's, that the dominant exchange in a high energy collision at fixed t is the 
pomeron. From what little we know about the dynamical make-up of the pomeron, 
we suspect it is some sort of collective excitation made of several gluons plus perhaps 
quarks. Thus to study the low energy regime of QCD, one way is to focus of particular 
bound states and try to derive their properties from QCD. However another option is 
to study the Regge families of hadrons, such as the pomeron, and try to calculate their 
parameters from QCD. We can not offer any reason why the latter option is better than 
the former. However, the one evident fact is that we have much better experimental data 
about Reggeons than about individual hadrons. Also from the point of view of light-cone 
kinematics, the description of Reggeons is more natural than of individual particles. If one 
accepts this line of reasoning to its fullest extent, one could imagine calculating masses 
of hadrons using Reggeon concepts. At present we do not have sufficient control on the 
approximations involved in our evolution equations to justify such calculations. However 
one could assume the radius of convergence for our equations is sufficiently large to make 
some sort of estimates. We will not discuss this point further in this talk. 

2. Evolution Equations 

The CLAP equation^, 

dni{x) , , sr^ f^ dx' 

where, 






e = -ln(l + -cln^) (2) 

c 7r /i^ 



and 



UJi 



V / dxuJi-,k{x), (3) 

^ Jo 



k ^0 



determines the Q^-evolution of the parton distributions ni{x), where x = ^ is the ratio 
of the parton to hadron longitudinal momentum in the light-cone frame. The splitting 
kernels, uji^k{-^), describe the inclusive probabilities of the parton decay into the opening 
phase space d^. Mellin transforming in In^ gives the anomalous dimension matrix 7(j) 
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for the twist two operators in QCD. For example in the case of pure gluondynamics, 

n,(x) = / ^(-)^e^^(^) (4) 

where 

2 1 

and ^{j) = Yn • Note that 7(j) = 'y{\—j)—'nctg{Tij) and 7(2) = due to the conservation 
of the stress tensor T^^. From eq. (4) we obtain that for a; ^ 0, ng{x) ~ ^exp{cJ ^\n. i). 

This imphes that total cross-section o"i(a;, Q^), for 7*p scattering at large energies y^ = -^ 
grows more rapidly than any power of Ins. This is a violation of the Froissart bound 
at < cln^ s. 

At small X, for parton distributions ng{x, kj_) depending on transverse parton momen- 
tum k±, one should use the BFKL equation^, 

dn{x, kj_) ^ 2uj{-kl)n{x, k^) + f d^k'^K{k^, k'^)n{x, k'^) (6) 



dln- 

X 

where n{x) = J d'^k±n{x,k±) and 

where k'j_ > 0. Note that the gluon Regge trajectory j{—k'j_) is related to u by j{—k^) = 
1 -|- tc;(— A;^). The kernel K for SU(A^c) gauge theory is, 

A-(.,..;) = |!/v.^^. (8) 

Observe that the infrared divergences cancel in the right hand side of eq.(6). 
The solution of eq. (6) can be written in the form^, 

n{x,k±) = - Y. e^'W du{-r'^'''"'hf'cm,u, (9) 

X ,y^ ,^- J — CO X 



m=—oo 



ki. 

Wh 



where Cm,u is determined by the initial conditions for n{x^ kj_) at fixed x, = arctg{ 
and the eigenvalue uo{v,m) of the corresponding stationary equation is, 

0^ fi dv -i+|m| 0^ 1 \in\ 

a;(z/, m) = ^N, J^ -^[y^ cos{u\ny) - 1] = ^iVc(^(l) - i?e^(- + ti^ + ^)).(10) 
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The biggest value of oo is C(j(0,0) = ^Nc\'ii'2, and therefore from eq. (9) we obtain that 
n{x,k±) ~ i(i)'^(0'0). This means that the solution of the BFKL-equation also does not 
agree with the Froissart bound. For this equation as well as for the GLAP-equation, 
the reason for this violation is that the evolution equations were obtained in the leading 
logarithmic approximation, where the S-matrix does not satisfy unitarity^. 

Thus we find in both cases, the GLAP and BFKL equations, the result is incomplete. 
As such we will construct a modified leading logarithmic approximation (LLA) that is 
compatible with the Froissart bound. 

3. Partonic Wave Functions 

The partonic distributions ni{x,kj_) are proportional to the imaginary part of the 
scattering amplitude at the momentum transfer q=0. It is natural to generalize the 
evolution equations for arbitrary momentum transfer. In this case the resulting equations 
could be considered as equations for the hadronic wave function. The usual Schrodinger 
equation Eip = Hip determines the mass of the hadron as a function of its spin, m^ = 
m'^{j). To determine j=j(^), one can replace this equation by the BFKL equation 

2HuiJ = Eij, (11) 

where 

and j = 1 + a; is the position of the j-plane singularity of the t-channel partial wave. 
The high energy asymptotics of scattering amplitudes are determined by the eigenvalues 
of equation (11) as A(s,t) ~ g^+^it) ^ 'pj-^g eigenvalues cu could in general also depend 
on t=— g^, but due to the conformal invariance of the BFKL equation^, in LLA this 
dependence is absent. The operator H12 on the left hand side of eq. (11) is^, 

^12 = I piJip 12 ^1*^2 In Ipul'PiP; + h.c. + ln(|PinP2n " 4^(1), (13) 

I "ill "2! 

where pi2 = Pi — P2, Pr = Xr + iyr, the momenta Pr = i-^, and h.c means the complex 
conjugated expression. 

To unitarize the results of the LLA, one must generalize eq.(ll) for compound states 
with an arbitrary number of gluons. Such a generalization was done in [1]. Here we 
discuss the BKP-equation for the large Nc case. Thus we consider the equation. 

Hip = Eip (14) 
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with E as given in eq. (12) and where the Hamiltonian H contains only interactions of 
neighboring particles, 

m 

H=J2Hr,r+l. (15) 

r=l 

The pair Hamiltonian i/r,r+i acts on the coordinates r and r+1 of the gluons as given by 
eq. (13). 

Note that there is also a generalization of the GLAP-equation (1) for matrix elements 
of quasipartonic operators of high twist^. This B'F'KL-equation is also simplified in the 
region of large Nc- It takes the form of eqs. (14) and (15) with the pair kernel describing 
the evolution of twist two operators. The eigenvalues of this equation are proportional 
to the anomalous dimensions of the quasi-partonic operators whose contributions are 
important in the small-x region. We will not discuss this equation any further below. 

4. The BKP-Equation in the Large Nc Limit 

From eqs. (13) and (14) one can derive the holomorphic separability of the Hamilto- 
nian, which is a central property for our present discussion. Thus we can write 

H = H + H*, (16) 

where H and H* act on the holomorphic {pj) and antiholomorphic (p*) coordinates re- 
spectively with 

n 

H = J:H,,+i, (17) 

i=i 

and similarly for H*. The pair holomorphic Hamiltonian is, 

Hj,j+i = ^ HPj,+i)Pj + -j^ ln(p,,,+i)P,+i + ln(P„ P,+0 - 2^(1). (18) 

An important outcome of holomorphic separability is that the solution of eq. (14) sepa- 
rates as , 

^(Pl, P2, ■ ■ ■ , Pn) = IZ V^(Pl, P2, ■ ■ ■ , Pn)V^(Pl, P2, " " " , Pn) (19) 

where the sum is over all degenerate solutions of the Schrodinger equation in the holo- 
morphic and antiholomorphic subspaces, 

E = e + i,H^ = eij, H*ij = i^j (20) 
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The pair Hamiltonian Hjj^i in eq. (18) can also be written in the forms®, 

H,,+, = HpI^.Pj) + ln(p^+iP,+i) - 2 lnp,,,+i - 2^(1) (21) 

= p,,,+iln(P,P,+i)p,"l+i + 21np,,,+i-2V^(l) (22) 

From the above representations, it is obvious that H is invariant under the Mobius 
transformations^, 

ap^ + b 

where a,b,c, and d are arbitrary complex parameters. The generators of these transfor- 
mations are 

n 

M = J2M„ M,' = pA M- =d^Mt = -p% (24) 

One can also obtain the transposed Hamiltonian H^ from H by two different similarity 
t r ansf ormat ions , 

H'^ = P,P^...P^HP,-'P^'---P-' (25) 

= Pl2p2i---pnlHpi2p23---pnl- (26) 

This implies that there are two different normalization conditions for the solutions of eq. 
(14) which are compatible with eqs. (25) and (26). These are, 

\\H\'=Jrf[dprPri^ (27) 

r=l 

\\H\'= frfi^^i^. (28) 

r=l Pf^r+l 

From eqs. (25) and (26) we conclude that there is a nontrivial differential operator®, 

A = Pl2P2Z---PnlPlP2---Pn: (29) 

which commutes with H, 

[A, H] = {) (30) 

Below we will show that there are an infinite number of operators that commute with H. 

5. Equivalence Between Multi-color QCD at High Energies and an Exactly 
Solvable Spin Model 
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We can write down the operator A in eq.(29) as follows^, 

A = rtr(MiM2---M„), (31) 

where Mj is the 2*2 matrix constructed from the Mobius generators Mj in eq. (24), 

M. = ( P% ^\ \ (32) 

In representation (31) the operator A can be interpreted as a transfer matrix for a lattice 
theory. On the hnks in the "space" direction (the auxihary subspace) there are discrete 
variables ^ taking values ^ = ±1 and on the links in the "time" direction (the quantum 
subspace), there are continuous variables p. 

To solve eq (20) exactly, one should find the one parameter family of integrals of 
motion, including the operator A of eq.(31). It turns out^ that such a family is the 
following, 

t{e)=tr{L^{e)L2{9)---Ln{e)), (33) 

where, 

is the so called L-operator. Let us also introduce the monodromy matrix, 

T{e) = Li{e)L2{e)---Ln{e). (35) 

One can verify^ that it satisfies the following Yang-Baxter equation: 

T^H' {u)T^'^' {v){u -v-Pi2) = {u-v + Pu)T^'^'{v)T''^'{u), (36) 

where Pu is the operator that interchanges the matrix spin indices (the right and left 
ones correspondingly). By taking the traces over indices ir and ir', we obtain: 

t{u)t{v) = t{v)t{u), (37) 

so that the operators defined in eq. (33) commute with each other. 

Now we want to prove that the operator t{6) of eq. (33) also commute with the 
holomorphic Hamiltonian in eq. (20). For this purpose the idea we use is^ that the spin 
model with the transfer matrix (33) can be considered as a modification of the Heisenberg 
model. However, instead of the fundamental representation of the group SU(2) with spin 
S=|, here we have the infinite-dimensional representation of the Moebious group SU(1,1) 
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with spin S=0 . For this new spin model, there is an unique Hamiltonian describing the 
interaction of nearest spins for which the model is exactly solvable. The general method 
of obtaining this Hamiltonian was developed many years ago^. Briefly, to do this one 
should construct the operator L^'^{6), which satisfies the trilinear Yang-Baxter equation 
for the case when both the quantum and auxiliary subspaces are one-dimensional (pi and 
P2). Then for this new spin model, the Hamiltonian is given by eq. (17), where Hi 2 can 
be calculated from the small-^ expansion of L^'^(^): 

L''\9) = P'^^il + 9Hi,2 + ---). (38) 

Here P^'^ is the operator which interchanges the coordinates pi and p2- According to the 
general theory^, L^''^{6) should also satisfy the linear Yang-Baxter equation: 

Li{u)L2{v)L^'\v -u) = L^'\v - u)L2{v)Li{u) (39) 

In this equation, the Li operators are 2*2 matrices (34). From eq. (39) we find^ that Hi^2 
depends only on the Casimir operator of the conformal group. This can be written in the 
form: 

(Ml + Ms)^ = m(m - 1) (40) 

We also find that i/1,2 satisfies the equation 

[i7i2(m)-/Ji2(m-l)](m-l) = 2, (41) 

for which the solution is 

Hi^2=i^{m)+^{l-m)-2^{l) (42) 

up to an additional term A(m), which is a periodic function (ie. A(m) = A{rh + 1)). 
Using eq. (10), we can verify that the expression for Hi 2 determined by eqs. (18) and (42) 
coincide. Thus, according to the general theory in [8], the Hamiltonian (17) commutes 
with all operators of the type t{6) in eq.(33). 

6. Conclusion 

We have shown above that in the generalized leading logarithmic approximation, the 
equation for the compound states of n-reggeized gluons is significantly simplified in the 
large A'"c-limit. In particular, it is conformally invariant and the Hamiltonian has the re- 
markable property of holomorphic separability. In addition, the equations for holomorphic 
and antiholomorphic wave functions have a sufficient number of conservation laws to be 
exactly solvable. This is related with the fact that the Hamiltonians in the corresponding 
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subspaces coincide with the local Hamiltonians of the exactly solvable Heisenberg model 
for spin S=0. As such, the quantum mechanical problem posed in eq. (20) is reduced to 
the pure algebraic one of constructing the representations of the Yang-Baxter algebra in 
eq. (36). The simple method of finding these representations was developed in [10]. It is 
based on the solution of the Baxter equation 

A(A)g(A) = (A + 0"g(A + 2) + (A - iYQ{\ - i), (43) 

where n is the number of reggeized gluons, A(A) are the eigenvalues of the operator t{iX) 
in eq. (33) and the function Q is the integer function in the complex A- plane. The 
eigenvalue A(A) has the polynomial expansion in A, 

A(A) = 2A" - 3{3 + l)A"-2 + ■ ■ ■ + A, (44) 

where n is the number of reggeized gluons, m=j-l is the conformal weight of the corre- 
sponding composite operator, j(j+l) is the eigenvalue of the Casimir operator (I]jMj)^, 
and A is the eigenvalue of the integral of motion A. The eigenvalues and eigenfunctions 
of eqs. (20) can be expressed through Q{X). For n=2 eq. (43) is solved in terms of 
hypergeometric functions. For n=3 the solution of eq. (43) for integer j can be expanded 

(2) 

as a linear combination of its solutions Q} for n=2 as, 

Q{\) = J2d,{A)Qt\X), (45) 

fc=i 

where the parameter A is determined in eq.(44) and dk{A) are orthogonal polynomials 
satisfying the recurrence relations, 

Ad,{A) = ^^^^^(A; - j){k + j + l)(4+i(A) + 4-i(A)). (46) 

2k + I 

The quantization condition for the eigenvalues A is, 

d,{A) = 0. (47) 

It is possible to express the energies e in eq. (20) directly in terms of dk{A) , when eq.(46) 
is analytically continued to complex j. The solution of eq. (46) would give a possibility 
to find the Odderon intercept^^. 
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